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Abstract 

• CP violations in the kaon system are studied in supersymmetric SU(5) model with 

right-handed neutrinos. We pay a special attention to the renormalization group effect 
Qh! on the off-diagonal elements of the squark mass matrices. In particular, if the Yukawa 

couplings and mixings in the neutrino sector are sizable, off-diagonal elements of the 
right-handed down-type squark mass matrix are generated, which affect CP and flavor 
violations in decay processes of the kaon. We calculate supersymmetric contributions 
^ ' to e (as well as Amx), Br{KL t:^vi>), and e'/e in this framework. We will see that 

' the supersymmetric contribution to the e parameter can be as large as (and in some 

case, larger than) the experimentally measured value. We also discuss its implication 
to future tests of the unitarity triangle of the Kobayashi-Maskawa matrix. 



1 Introduction 



One of the most important issues in particle physics today is to understand the origin of 
CP violations. Indeed, many efforts have been made to measure CP violations in various 
processes. So far, for the K system, non-vanishing values of the e and e' parameters have 
been observed. In addition, CP violation in the ipK^ process, i.e., the angle 



IS 



now being measured by the on-going S-factories, and even at the present stage non-vanishing 
CP violation in this process is reported fl], ||. 

In the framework of standard model (SM), the most well-known mechanism to explain 
these CP violations is to introduce the Kobayashi-Maskawa (KM) matrix which con- 
tains one physical complex phase for the three family case. That is, using the Wolfenstein 
parameterization [ffl: 



VkM ^ 
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the parameter rj parameterizes the size of the CP violation. 

Importantly, the measurements of e and (pi are used to constrain the parameters in the 
KM matrix. Currently, all the measurements of the CP violations more or less suggest the 
same region on the p vs. t] plane, and hence the observed CP violations are well explained 
in the framework proposed by Kobayashi and Maskawa. In the near future, some of the 
measurements of the CP violations will become more precise ^]. In addition, there will be 
other constraints on the p vs. rj plane from various new processes like n^uu and decay 

processes of 3^ and Bg mesons. These processes will provide important tests of the unitariry 
of the KM matrix and the origin of the CP violations. In particular, if some new physics 
exists, it may provide a new source of the CP violation which may be seen as a deviation 
from the SM prediction on the p vs. 1] plane. 

Of course, possible deviations depend on featrures of new physics beyond the standard 
model. Among various models, in this paper, we consider one of the most well-motivated 
ones, that is, supersymmetric (SUSY) unified model with right-handed neutrinos. Such a 
model can solve some of the theoretical and experimental problems which cannot be solved in 
the framework of SM. First, in the supersymmetric theories, the serious naturalness problem 
can be avoided because of the cancellation of the quadratic divergences between bosonic 
and fermionic loops. In addition, in the minimal SUSY SM (MSSM), successful gauge 
coupling unification can be realized contrary to the standard-model case where three gauge 
couplings do not meet at any high energy scale. Furthermore, in this framework, solar and 
atmospheric neutrino problems 0, ^ may be solved by small neutrino masses generated 
by the seesaw mechanism |p. Thus, supersymmetric unified theories with right-handed 
neutrinos are theoretically and experimentally well-motivated, and it is worth studying its 



"^""^The angle is also called the angle (3. In this paper, however, we do not use this notation since we use 
the angle (3 to parameterize the vacuum expectation values of two Higgs bosons. 
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phenomenological consequences. In this context there have been some phenomenological 
studies fl^, |Tl| which obatained interesting results in some processes. 

In this paper, we study the CP violating processes in the kaon system in the framework 
of supersymmetrics grand unified theories (GUTs) with right-handed neutrinos.0 Indeed, 
in such models, there are many possible new physical phases in the soft SUSY breaking 
parameters, which may affect various CP violating processes. 

This paper is orgenized as follows. In Section ^ we introduce the model we consider. 
Then, in Section ^, rates of the various CP violations in the /^-system is numerically evalu- 
ated. Sectionals devoted for the conslusions and discussion. Relevant formulae are collected 
in the Appendices. 



2 Model 

In this paper we consider the minimal SU(5) GUT with singlet right-handed neutrinos. Let 
us denote 10, 5, and the right-handed neutino chiral multiplets in i-th generation as \E'j, $j 
and Ni, respectively. The superpotential is given as follows;^ 

Wgvt = l^i[Yu]ij^,H + ^i[Yn]^,<!>jH + N,[Y^],,<l>jH + In,[Mm],^N„ (2.1) 

o Z 

where H and H are 5 and 5 representation Higgs fields, respectively. Here i,j are generation 
indices. M^r is Majorana mass matrix of the right-handed neutrinos. Throughout the paper 
we adopt the universal structure of Mat for simplicity, 

[M^h = M,J,,. (2.2) 

For our discussion it is convenient to choose basis of \E'j, $j and Ni such that the down- 
type Yukawa matrix is diagonal and mixing matrices appear in the up-type and neutrino 
Yukawa couplings. The mixing matrix of the up-type Yukawa coupling corresponds to the 
KM matrix Vkm- In general the mixing matrix of the neutrino Yukawa coupling is given by 
a combination of two unitary matrices that are mixing matrices of left-handed neutrinos and 
right-handed Majorana neutrinos. In our case the mixing matrix of the neutrino Yukawa 



coupling becomes the Maki-Nakagawa-Sakata matrix I^mns |16| because of the universal 



structure given in Eq. ( |2.2| ). Hence in this basis Yukawa matrices are decomposed as 

Yu{Mgvt) = V^u^qYuVkm, Id(Mgut) = Y^, Y^{Mgvt) = Y^V^^s^l, (2.3) 
where F's are diagonal matrices; 

Yu = diag{yui,yu2, Vus), Yd = diag(?/di, l/d2, Vds), Yn = diag(y„i, yn2, Vns), (2-4) 

^^Flavor violation in K and B sysytenis in the framework of SUSY SM without the right-handed neutrino 
have been discussed in |l|, . 

"^^In addition, an adjoint Higgs S which is responsible for SU(5) breaking is assumed. However interactions 
between E and the other chiral multiplets are not significant for our discussion. Hence we do not discuss 
these interactions in this paper. 
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and G's are diagonal phase matrices; 

• (Q) • (Q) •(<?)- • (L) ■ (i) . (L) 

eQ = diag(e*^i ,e*^2 ^ g^'^^ ), 6^ = diag(e*'^i , e*^^ , g^v'a ). (2.5) 

Since physics does not change under redefinition of overall phases of Qq and G^,, we fix 
(f^^^ + (f^^^ + (fl^^ = and + ¥^2^^ + V's^^ = 0. Vkm and Vmns are parameterized by four 



parameters, i.e., three mixing angles and one phase as follows [|^; 



KM 
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(2.6) 



where c^-^^ = cos 6'^^'' and s^-^^ = sin6'|^'', and Vmns is obtained by exchanging Q ^ L. 

Let us decompose VFgut by using the standard model fields. We embed the standard 
model fields into the GUT fields as 

where (5i(3, 2)i/6, [7^(3, l)_2/3, A(3, 1)1/3, Li{l, 2)^1/2 and Ei{l,l)i are quark and lepton 
fields in i-th generation with SU(3)^ x SU(2)^ x U(l)y quantum numbers as indicated. Here 
a, P, 7 represent SU(3)(^ indices and a, b are SU(2)^ indices. With this embedding, W^gut at 
the GUT scale becomes 



W^GUT = WssM-lQ^[V;[M^QYuVKMhQJHc-E,[<^LYu],,U,Hc + uSl^^^ 

-Q^lYnQLhL^Hc + NilY^V^^^QLj^.D^Hc, (2.8) 
where He and He are colored Higgs fields, and low energy superpotential is given as 



WssM = H^ilYuVKuhQj + H^DilYnhQj + HdEi[YE\ijL 
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+N,[YnV^^^],,L,H^ + ^M.^iV^iV,. (2.9) 

Notice that the phase matrices 6q,l <io not appear in Eq.(p.9D. Hence these phases v^l*^'^^ 
are independent of the CP phase in the KM matrix. As we will show later, GUT phases 
(/?^^^ can have important implications to CP violations in the kaon system. Although simple 
SU(5) GUT predicts Yd = Ye at the GUT scale, this relation is unrealistic for the first 
and second generations. In order to explain realistic fermion mass pattern, some new fiavor 
physics are necessary. Although these new physics can provide extra source of fiavor and 
CP violations [^, |2^ , we do not consider these effects in this paper. 

With the superpotential given in Eq. ( |2.9| ), the left-handed neutrino mass matrix is gen- 
erated by the seesaw mechanism, 

[mvj]ij = ^j^/^ [^MNS^7V^MNs]ii = [^MNS^TV^MNslij' (2-10) 



where v ~ 246 GeV and tan (3 is the ratio of two Higgs vacuum expectation values. In 
our analysis, we use the masses and mixing angles of neutrinos suggested by solar and 
atmospheric neutrino data. Atmospheric neutrino data [0 at Super-Kamiokande implies 
that z/^ — mixing is almost maximal, hence we set s^? ~ 1/V2. As for Ue — mixmg, 
we consider large and small angle MSW solutions to the solar neutrino problem [21].^ For 
large angle MSW solution, we use 



/ 0.91 -0.42 < 1 

-0.30 0.64 0.71 I , (2.11) 
V 0.30 -0.64 0.70 



~ (0, 0.004eV, 0.03eV), 1/mns 
and for small angle MSW solution, 

~ (0, 0.003eV, 0.03eV), Hins ^ I -0.028 0.70 0.71 I . (2.12) 




[Vmns]i3 is small as indicated by CHOOZ experiment [^. In our analysis we set [Vmns]i3 = 
unless otherwise stated. Nonzero value of [Vmns]i3 can affect our result, especially for the 
small angle MSW case as we will discuss later. 

In our basis, Yukawa matrices for down-type quarks and charged leptons are diagonal at 
the GUT scale. Their diagonalities do not hold at the weak scale because of renormalization 
group (RG) effects. However these effects are minor for our discussion and Di and Lj fields 
in Eq. (|2.9| ) are mass eigenstates with good accuracy. Hence it is very useful to use this basis 
to see the flavor violations in K and B meson system caused by the right-handed neutrinos. 

Although mass matrices of down-type quarks and charged leptons are diagonal in our 
basis, down-type squark mass matrix may have flavor violating off-diagonal elements. The 
soft SUSY breaking terms for scalar fields are given by 

+ -'ii}i[Au]ij'ipjh + 'ii)i[AD\ij(i)jh + hi[AN]ij(i)jh, (2.13) 

o 

where ipi, (pi, -hi, h and h are scalar components of the superfields \E'j, $j, Ni, H and H, 
respectively. Soft SUSY breaking terms for the low energy effective theory below the GUT 
scale are given by 

^SSM = Qil'iTT'QhQj + ~di[m^D]ij~d] + "^iVnyiijU] + l\m]]ijlj + li[m|]jji] + n[m\]h] 

+Kui[Au]ijqj + hddi[Ad]ijqj + hdei[Ae]ijlj + Khi[An]ijlj, (2.14) 



There are other solutions to the solar neutrino problem, LOW, vacuum and Just-So solutions. In these 
solutions, second generation neutrino mass is smaller than that in large or small angle MSW solutions. 
m^^ ~ 3 X lO^^'eV, 1.2 x lO^^eV and 2 x lO^^eV for LOW, vacuum and Just-So solutions, respectively. 
When we impose perturbativity of y„3, j/n^ is very small because of the universality of Mjv- Therefore when 
s^3^ = 0, flavor violation for the kaon system is too small to be observed. So we do not discuss the three 
cases. But if is not universal or is nonzero, the situation may change. 
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where qi, Ui, di, k, Cj, hi hu and hd are scalar components of the superfields Qj, t/j, Z)j, 
Lj, Ei, Ni, Hu and if^^, respectively. With the embedding Eg. ( p.7|) the above scalar mass 
matrices at the GUT scale are given by 

= ^lo^^ = 0q'^km"^?o'^km0Q, ml = mf, 

ml = BLmlel ml = eLml^el (2.15) 

We are particularly interested in the RG effect on the off-diagonal elements of the 
sfermion mass matrices. Thus, in our analysis, we use the boundary condition of the mini- 
mal supergravity model. Then, the SUSY breaking parameters at the reduced Planck scale 

~ 2.4 X lO^^GeV are give by 



m 



Au{M,) = aoYuiM,), AoiM,) = aoYniM,), An{M,) = aoF7v(M,), 
Mi(M,) = M2(M,) = M3(M,) = Ml/2, (2.16) 

where Mi, M2 and M3 are SUSY breaking masses of U(l)y, SU(2)l and SU(3)c gauginos, 
respectively. Importantly, sfermion masses are universal at the reduced Planck scale. Below 
the Planck scale M^,, however, off-diagonal elements of scalar mass matrices are generated 
by the RG effect and sizable CP and flavor violations may be possible. Of course, if there 
are tree- level off-diagonal elements of the sfermion mass matrices at the scale M^, they may 
affect the rates of the CP and flavor violating processes. However, it is fairly unnatural 
to assume cancellations between the tree-level and RG contributions to the off-diagonal 
elements. Thus, we believe our approach will give us a conservative estimation of the rate 
of the CP and flavor violating processes unless some accidental cancellation happens. 

First we see lepton flavor violation. In the presence of the right-handed neutrinos, Hu — 
N — L interaction, the fourth term in Eq.( p.9|) , violates lepton flavor. Although the right- 
handed neutrinos are very heavy, this interaction affects the RG evolution of slepton mass 
matrix and generates off-diagonal elements of m| at the weak scale. Such an effect can be 



probed by ;U — > 67, r — > /i7, /i — e conversion, and so on 123, 121, 



Since down-type quarks and charged leptons belong to the same representation of SU(5), 
flavor of down-type quarks is also violated via the right-handed neutrino Yukawa interaction 



lOl , |TT[]. Indeed interaction among N, D and He (the last term of Eq.( p.8| )) induces off- 
diagonal elements of the right-handed down-type squark mass matrix, which are evaluated 
with the one iteration approximation as 

[mlh ^ ^{Sml + a^)e^(^-''-^^'') E?/4[VMNs].4V^MNs].;^log ^ i^ ^ j)- (2.17) 

Notice that the size of the off-diagonal elements is determined by the parameters in the 
neutrino sector, i.e., neutrino Yukawa coupling and Vmns- In addition these elements depend 



"^^In the GUT framework, lepton flavor is violated by interaction with colored Higgs. Such a elTect have 
been studied in il9|, M EtII 
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on phase (yj,-^^ — V'j^'*. Hence flavor violation and new phases caused by the right-handed 
neutrino affect flavor and CP violations in the quark sector. In this paper we focus on kaon 
system to investigate such effects. 

In order to discuss flavor violation, it is useful to introduce the following variable which 
is off-diagonal elements of m|i normalized by squark mass scale rriq, 

Alf ^ l!% = t^. (2.18) 

Here we choose the first generation down-type squark mass as squark mass scale, which is 
approximately estimated as ~ ml + Y.IM^^ our model. Using S23'' — l/V^ and 

(■^13^)^ 1) ^12^^ is approximately estimated as 

A[f ~ 5.5 X 10"^ 

Uoi^CevJ [omeVj + W, '' J ^ 

From the equation, we find that the flavor violation is enhanced for larger M^,^. It reflects one 



characteristic feature of the seesaw mechanism, i/n^ oc M^,^ with fixed left-handed neutrino 
mass. 

A^f^ is sensitive to the neutrino mixing angles. First we discuss s[^^ = case. For the 
large angle MSW solution, A^f'' is not suppressed by mixing angle because 5^2'* ^12'' ~ 1/2- 
On the other hand, Aj^f-* is suppressed by the factor Si2^Ci2 ~ 0.04 for the small angle 
MSW case. Hence, for fixed M^,^, amphtudes for flavor violating processes in the large angle 
MSW case are larger than those in the small angle case by one order of magnitude. Next 
we consider an effect of finite s[^\ Since mj^.^/mj^^ ~ 10, the effect becomes important when 
~ 0.04 for the large angle MSW case, and s^^^ ^ 0.004 for the small angle MSW case. 
Therefore improvement of the upper bound on 3^3'' can have impact on the quark sector 
flavor and CP violations, especially for the small angle MSW case, in our model. 

Off-diagonal elements of the left-handed squark mass matrix are generated by quark 
flavor mixing via I^m, 

\<h = ^(3mg + al)y^,[V^y,UV^u\,, (s log ^ + log ^) ^ j). (2.20) 

As in the dn sector, we introduce the following variable 



^L) ^ (2.21) 



A^2'' is approximately given by 



Ai^) ^ 1.6 X 10-4 1 3mg + ag ^,^^^ 



sm p rriq 



Comparing the above expression to Eq.(p.l9D, we find that dominates over A^2 when 



M,^ > 10^4 GeV. 
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3 Numerical Results 



3.1 e and Amx 

In this subsection we discuss effects of the right-handed neutrinos on the AS* = 2 processes. 
The effective Hamiltonian for the AS* = 2 processes is given by 



1=1 i=4 

where operators are 

Ql,1 = 4:{do,J^,PLSa){dl3YPLSp), 

Ql,2 = A{daPLSa){dpPLSp), 

Ql,3 = 4:{daPLS(3)idf3PLSa), 

= 4:{daPLSa){dpPRSl3), 

Qs = 4((i„PiS^)(J^PRS„), 



(3.1) 



(3.2) 



and the operators (z = 1 — 3) are obtained from Qi^i by exchanging L R. Here a and 
(3 are color indices. In order to calculate Amx and e, we take account of QCD correction 
to the Wilson coefficients, and using formulae for the leading-order QCD correction given in 
the matrix element of the effective Hamiltonian at /ic = 1-3 GeV is given as 



+ 
+V5C5 



rriK 



1 1 

12 2 Vm^ -|- rrid 



1 + i ( 

4 6 Knis + rriii 



(3.3) 



where we used vacuum saturation approximation to calculate matrix elements, and the QCD 
correction factors are: 



a,(mc) 



6/25 



6/23 



0.77, 



(3.4) 



and ri4 = rii ~ 2.8 and 775 = rji ~ 0.88. 

With this matrix element, e and Attik are given by 

ei^/Hm{K^\nes\K'^) 
2\/2mKAmK 

AniK = —\{K^\n,s\K^)\. 
rriK 



(3.5) 
(3.6) 
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There are four kinds of the SUSY contributions to the effective Hamihonian, 1) gluino 
mediated, 2) chargino mediated, 3) gluino-neutrahno mediated and 4) neutrahno mediated. 
Among them, gluino mediated 0{al) contribution dominates over the others in our modeL 
Hence for a while we concentrate on the gluino contribution for an intuitive discussion. At 
the scale where SUSY particles are integrated out, Wilson coefficients are given in terms of 
mass insertion (MI) parameters as follows [0: 



m|, m|, m?, m?, m?, m?) 
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+ ll/g^(m|, m|, m?, m?, m?, m?)], (3.7) 



47r o/.l'P'lX^ A r . 0^0/ 2 2 2 2 2 2\ 

(m^, m^, mg, m^-, m^-, m^-J 



+ 11/6 (m|, "^1^ m?, m?)], (3.8) 

^4|gi = (Ai?) (A(?) m|[7m|/°(m|, m^-, m?, m^-, m|) 

7-1 / 2 2 2 2 2 2\i /o 

-/g (m^, m^, m^-, m^-, m^, m^)J, (3.9) 

~ Q IGvr^ 2 / A {L)\ ( A (_R)\ 4r2r0/2 2 2 2 2 2\ 

C^sImi = -^"s (^^12 j (^Ais'j m^-[m^/6("^G,"^G,"^g,"^g,"^g,mg-) 

+5/g (m|, m|, m?, m?, m?, m?)], (3.10) 

where the function is defined in the Appendix 0. Since flavor violating left-right mixing 
is very small, so are Cl,2\mi^ C*i?,2|Mi5 C*l^3|^j, Cr,z\mi- Therefore we neglect them in the 
following rough estimation. Although these coefficients are suppressed by l/m.^, they can 
be comparable to the SM contributions since they are 0{a1) while the SM contributions are 
0(q;|). Notice that contributions from operators Q_r,i, Q4 and Q5 are sizable only when 
flavor in sector is violated. 

The flavor and CP violations in sector have important implication to Am^ and e 
because the matrix elements of the LR type operators Q4 and Q5 are enhanced by the factor 
{mK/TTisY ~ 10. In addition the Wilson coefficient of the operator Q4 at /ic is enhanced 
by the QCD correction factor 774. Hence contribution from the operator Q4 dominates over 
the other SUSY contributions when and A^f'' are comparable. This is the case for 
the minimal SU(5) GUT with the right-handed neutrinos, and we shall see that SUSY 
contribution to e can be as large as the experimental value. In the absence of the right- 
handed neutrinos, the LR type operators do not exit, and SUSY contribution to e and Amx 
is more suppressed. 

Before showing numerical results, let us discuss how e and ArriK depend on model pa- 
rameters. To obtain approximate formula, we consider only the contribution from Q4 since 
this operator has the most important effect. In addition, for simplicity, we consider the 
case where the gluino is as heavy as squarks. Then SUSY contributions to e and Am^ are 
estimated as 



2 

V nig ) VlOi^GeV; V0.004eV; 
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x"°('^'":ff'+^'). (3,11) 



sin P 

^ \ J VlQi^GeVy V0.004eVy 



x-^f4?cS^^ + ^4?). (3.12) 
sm p y J 

From the above estimation, we find that e^^SY large as experimental value e'^^^ = 

(2.271 ± 0.017) X 10-3 if M^^ > lO^^GeV, for ruq ~ ITeV. On the other hand Amf^^ 



is much smaller than the experimental value Am^^ = (3.489 ± 0.008) x 10-^^MeV [|T] 
unless M,y^ is larger than about 10^^ GeV. We do not consider such a large right-handed 
neutrino scale since the neutrino Yukawa coupling [lArJss blows up below the Planck scale for 
M^^ > 2 X lO^^GeV. Hence the CP violat ion in dji sector can affect e and does not contradict 
to ArriK measurement. 

Now we show our numerical results. In the numerical calculation not only the gluino 
contribution, but also other contributions with chargino and neutralino propagators are 
included. Furthermore contributions from all eight operators shown in Eq. (|3.1|) are taken 
into account. Numerically we checked that the gluino contribution through the operator Q4 
is dominant SUSY contribution. 

In Fig. |I|, we show e^usY ^(l) _ ^(l) _^ plane for large angle MSW solution 

with M^^ = 3 X lO^^GeV, tan/5 = 3, M2 = ISOGeV and = 0. From the figure, we 
find that e^^^"^ is maximized at 93^^^ — + 0i ~ 7r/2, as suggested by the approximation 
Eq.( |3.1lD . Notice that heavy Majorana neutrino mass and large mixing angle induce large 



lepton flavor violation in e and /i sector as well in our model. Since such lepton flavor 
violation is most severely constrained by ;U 67, we also calculate Br{^ 67) assuming 
Yd = Ye at the GUT scale as the simple SU(5) model predicts. For the above parameters, 
Br(fj, 67) is larger than the experimental upper limit 1.2 x 10"^^ unless rriq ^ 980GeV as 
shown in Fig. |l|. g^^SY ^^^Qgg ^q^^ depend on tan (3 so much. On the other hand, the constraint 
gets severer for larger tan f] since the amplitude for /i — 67 is almost proportional to tan /3. 
Notice that, even with relatively heavy squark mass of irig ^ 980 GeV which is consistent 
with the fi ^ e'-f constraint, the SUSY contribution may be as large as (or even larger than) 
the experimentally measured value. In addition, the fi —>■ e'-f constraint here is based on 
the relation Yd = Ye at the GUT scale, this relation may be violated by some new flavor 
physics which gives realistic fermion mass pattern for the first and second generations. Such 
a mechanism can change the prediction for /i ^ 67 branching ratio from that in the simple 
SU(5) case [|1^, |1^, Hence we should keep in mind that this constraint have uncertainty 



stemming from nontrivial texture of the fermion mass matrices. 

In Fig. we show the dependence of e^^SY squark mass rrig in the large angle MSW 
case, for M^^ = 5xl0^^ 1x10^^, 2x10^^ and3xl0i^ GeV. As one can see from Eq. (gATJ) , the 



higher the right-handed neutrino mass scale is, the larger e^usY bg^omes. If M,^^ ^ 10^^ GeV, 
gSUSY jg comparable to the experimental value as estimated before. It is easy to understand 
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Figure 1: Contour of le^usYj units of 10 ^) on vs. {(pi V^2 0i)/^ plane for the 
large angle MSW case (solid lines). We take M2 = ISOGeV, M^^ = 3 x lO^^GeV, Oq = 
and tan (3 = 3. In the left side of the dashed line, Br{fi 67) is larger than experimental 
limit 1.2 X 10-", provided that Yd = Ye at the GUT scale. 



the dependence on rriq. When mo is large, A^^^ is also enhanced. On the other hand when 
squarks become much heavier, they decouple from low energy physics and e^usY ^^^^ smaller. 
Hence, for M2 = 150 GeV corresponding to M3 ~ 450 GeV, e^usY ^g^j^gg ^]^q maximum value 
at rriq ~ 800 GeV and monotonously decreases for rriq ^ 800 GeV. 

The same dependence in small angle MSW solution is shown in Fig ^, for M^^ = 1 x 10^^, 
2 X 10^^ and 3 x lO^'* GeV. Solid lines correspond to 5^3'' = case. In this case e^^^"^ is 
much smaller than that in large angle solution since S12 is very small. However if Si^^ is 
nonzero, the situation changes drastically. Dashed lines show results with 5^3'' = 0.02. In 
this case g^usY jg enhanced and can be comparable to e^^^ since rriiy^ is larger than rrii^^ by 
one order of magnitude and 532'' is also large. As 5^3^ gets close to the current upper limit of 
about 0.15, e^usY i^g^omes much larger. ^ Hence improvement of the limit on 5^3'' obtained 
by reactor experiments is very important to discuss CP violation in the kaon system in our 
model, especially for small angle case. 

It is important to consider the implication of the SUSY contribution to e to the KM 
matrix since in the standard model, the e parameter provides an important constraint on 



the p vs. T] plane. In the standard model, the e parameter is given as (see, for example, |^ 



e^^ ~ C^BkA^X't] -ViSoix,) + VsSoix,, Xt) + A'X\l - p)v2Soixt) e'^'^ (3.13) 



The effect of finite sjg^ on /x ^ 67 has been discussed in 
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Figure 2: e^^^^ in tlie large angle MSW case for M^^ = 5 x 10", 1 x 10^^ 2 x 10^^ and 
3 X lO^^'GeV from the below. We take M2 = 150 GeV, oq = 0, tan /3 = 3. The horizontal 
dash-dotted line corresponds to experimental value — 2.271 x 10~^. 
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Figure 3: e^^SY ^^^^ ^^g^g ^^^^^ f^j. ^g) ^ q ^gQ^-^ ^-^g-j ^g) ^ q q2 (dashed 

hne). In both cases, we choose M^^ = 1 x 10^"^, 2 x 10^^ and 3 x 10^"^ GeV from the below. 
Wc take M2 = 150 GeV, ao = and tan P ~ 3. The horizontal dash-dotted line corresponds 
to experimental value e^^P = 2.271 x 10~^. 
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where = 3.78 x 10*^, t]i — rj^ are parameters for the QCD factors given by rji = 1.38 ±0.20, 
r]2 = 0.57 ± 0.01, and rj^ = 0.47 ± 0.04, and Inami-Lim functions |3^ So{xc,t) and So{xc,Xt) 
are given by 



So{Xc,Xt) 



4:Xc,t -llxl^ + r 



c,t 



4(l-a;,,t)2 



3x1^ \ogXc,t 
2(1 -s. 



c,tl 



Xn 



3xt 



?>x1 log Xt 



4(1 -Xi) 4(1-Xi)2 



(3.14) 
(3.15) 



with Xc,t = ^Itl'^^w- addition, we use Bk = 0.75 — 1.10 |T^. Comparing Eq. ( p.l3| ) with 
the experimentally measured value of e, a constraint on the p vs. r] plane is derived in the 
standard model. 

If SUSY contribution to e exists, such a constraint should be reconsidered. In particular, 
as we have seen, e^^^^ may be as large as the experimetally measured value of e, and hence 
the observed value of e may have a significant amount of a contamination of e^usY gj^j^g ^j^g 
size of e^usY jg model dependent and hence is unknown, we should disregard the constraint 
from the e parameter. Of course, the p vs. rj plane is still constrained by other quantities 
like I [VKM]Mfe/[^KM]cfe| and Am^^, which are related to the parameters in the KM matrix as 

|[VKM]«6/[VKM]cfe| ^ A(p2 + 



^1/2, and 



■r7ijmB^(F^ SijJm^S'o(xt)A^A2 



(3.16) 



where FbJBb^ = (200 ± 40) MeV, and r]B = 0.55 ± 0.01 Ip 



Importantly, the allowed region on the p vs. rj plane changes as we exclude the constraint 
from e. Accordingly, possible size of the CP violation in the KM matrix, i.e., the //-parameter, 
changes. 

To make a more quiantitative discussion, we derived constraints on the p vs. 1] plane 
comparing the standard model predictions on e, | [Vkm]u6/[Vkm]c6| and ArjiB^ with their 



experimental values.^ With and without the constraint from e, we found the upper and 
lower bounds on the rj parameter become^ 



0.26 <ri< 0.48 
0.17 <ri < 0.53 



with e, 
without e. 



(3.17) 
(3.18) 



where we assumed flat distributions of the uncertainties in various parameters given above, 
and we adopted A = 0.834 ± 0.039 [|17|. As one can see, the upper and lower bounds on 
7] changes by about 10 % and 35 %, respectively, if the constraint from e is not taken into 
account. This fact has important implications to future studies of CP violations using rare 
processes like — > tt^uu, as will be discussed in the following subsection. 

"^^Upper bound on Ams^ also provides another constraint. Inclusion of such information, however, does 
not change the following result significantly, and hence we neglect it in our analysis. 

'^^For the case without e, negative value of 77 is also allowed. However, < is strongly disfavored by the 
BELLE and BABAR experiments |l, and we do not consider such a case. 
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3.2 Kl TT^Z/Z^ 

In the future, we may have another interesting information of the CP violation via the rare 
process Kl tt^uu. In the standard model, this process is mostly a signal of direct CP 
violation, and the amplitude for this mode is proportional to rj. Most importantly, hadronic 
contribution to — > n^uu is small, and short distance effects of QCD are well under 



control [Q. Therefore, this mode is theoretically very clean. Although the branching ratio 
for this process is very small, we may have substantial number of —>■ tt^uu event in future 
experiments [Q, which may provide important test of the unitarity of the KM matrix. Due 
to Ref. 0, the rj parameter may be determined at 10 % accuracy. 

This fact has an important implication to the search for the new physics. In the standard 
model, the t] parameter determined by — > n^uu should be consistent with that from other 
constraints, in particular, that from e. If the SUSY contribution to e is sizable, however, this 
may not be realized. Indeed, as discussed in the previous subsection, the upper and lower 
bounds on rj may vary by about 10 % and 35 %, respectively, if the constraint from e is not 
taken into account. Given the fact that r/ may be determined with the accuracy of 10 % 0, 
measurement of rj in the experiment of Kl — >■ ir'^uu gives an important constraint on p vs. 1] 
plane independetly with the measurement of e. 

Of course, the Wilson coefficients for Kl tt^uu may be directly affected by the SUSY 
loops. Thus, we also calculated the SUSY correction to the Kl — * 7r° z/z/ process in our 
framework. 

There are penguin and box diagrams in SUSY loop corrections which induce Kl — > 
TT^uu. Effective operators are Ql^ = {dLYsL){i'Llu^L) and Q^^ = {dRn^SR){vLlvi^L) with 
corresponding Wilson coefficients Cll and C^^. SUSY contributions to Wilson coefficients 
are summarized in the Appendix p. For the process Kl tt^uu, we found that the total 
SUSY contribution to the Wilson coefficients is a few % of the SM contribution, and that the 
SUSY contribution has distractive interference with the SM one in our model. As a result, 
it is challenging to see this deviation in the experiments. 

For reasonable choices of parameters, the dominant SUSY contributions are from chargino 
and charged higgs penguin diagrams. Thus, the CP violating phase is controled by the phase 
in the KM matrix since the SUSY contribution is approximatelhy proportional to the 1-2 
element of the left-handed squark mass matrix. On the contrary, effects of the GUT phases 
are tiny in Kl — >■ vr^i/z/ because of the smalless of neutralino contribution which is, we 
found, mostly 0.1 % level. Since the SUSY contributions to this mode are quite small to 
be observed, useful information about the parameter rj is expected from Bt^Kl —>■ tt^uu) in 
our model. 



3.3 



e e 



In supersymmetric models, the parameter e' can be also modified due to the supersymmetric 
loop effects. In particular, in Ref. [0], it was pointed out that the SUSY contribution to e'/e 
may become much larger than the experimental value, [e'/eY^^ = (19.3 ± 2.4) x 10~^ [^ , 
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because of the significant modification of the A/ = | amphtude contributing to e'/ e.|^ As a 
result, in some parameter space, e'/e provides severer constraint on CP and flavor violations 
in the soft SUSY breaking parameters than the e does, and hence one might worry if the 
SUSY contribution to e'/e is in a reasonable range in our framework. In this subsection, we 
discuss the SUSY contribution to e'/ e. 

Before discussing the supersymmetric effect on e'/e, let us first introduce several formulae 
necessary to evaluate e'/e. The e'/e parameter is calculated once the Wilson coefficients for 
the AS* = 1 effective Hamiltonian are given. The AS* = 1 effective Hamiltonian is denoted 



as 



where Cj are Wilson coefficients, and 



^(AS- 


=1) 


= A{sa-1^,PLdo)Y^{qpl^,PLqp), 
<? 


(3.20) 




=1) 


= A{sali,PLdp)Y^{qp-f^,PLqa), 
<? 


(3.21) 




=1) 


= A{sal^,PLda)Y.(^P^^^PR(lp)^ 
1 


(3.22) 


Qtt 


=1) 


= 4:{sc,-ff,PLdf3)J2iQl^^i^PR(la), 


(3.23) 




=1) 


= 4:{sal^PLda) E (^qi^pltiPaqp), 
g 


(3.24) 




=1) 


= 4:{sal^^PLdf3)J2eq{qf3-ff,PRqf3), 
<? 


(3.25) 


^{AS- 


=1) 


= 4:{sc,-f^PLda) J2 eq{qp7t^PLqp), 


(3.26) 




=1) 


= A{sa-i^,PLdp)^eq{qp'-i^PLqo), 

9 


(3.27) 



with Pl/r = |(1=F75)5 « and (3 being the color indices, Cq denoting the quark electric charge, 
and Q^^i~^^ = QL^f~^'^\L^R- Based on the structure of the operators, we call Q^^f~^^ with 

z = 3, 4, 9 and 10 {Q^r^~^^ with z = 3, 4, 9 and 10) as LL-type (i?-R-type) operators while 
others as LR- or i?L-type operators. 

In the standard model, the gluon-penguin diagram generates the operators Q^^f~^^ with 
i = 3 — 6, and hence their Wilson coefficients are of O(awas) where aw indicates the 
coupling constants for the electroweak interaction. On the contrary, the operators Q^^f~^^ 

^^In SUSY models, e'/e can be also modified if the A-parametcrs for the down-type quarks are not aligned 
to the down-type Yukawa matrix p6| . 

#iOWe use the suffix {AS = 1) for the Wilson coefRcients and operators for the AS* — 1 effective Hamiltonian 
to distinguish them from those for the AS = 2 effective Hamiltonian. 
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with i = 7 — 10 are only from the electroweak processes and the corresponding Wilson 
coefficients are of O(a^). Although the Wilson coefficients contributing to the AI = | 



amplitude are smaller than those contributing to the AI = | one, both amplitudes are 



significant since the AI = | contribution has extra enhancement factor. 



From the effective Hamiltonian given in Eq. ( p.l9D , e'/e is given by (see, for example, 
1) 

UJ 



e 2|e|ReAo 

+ (L R), (3.28) 
where |(7r7r)/) denotes the two-pion state with isospin /, and 

^ = ^5-^' 3.29 
Re^o 

with 

Ao = ((7r7r),=o|Kff|i^), ^2 = ((7r7r),=2|Kff|^). (3.30) 

Numerically, 00 ~ 0.045. In addition, Q^jri' is from the isospin breaking in the quark masses, 
and numerically we use = 0.25 ||31|. Since is a small number, the AI = | contribution 
in Eq. ( p.28|) , which is proportional to 00^^, is enhanced relative to the A J = ^ contribution. 

Now, let us consider the supersymmetric contribution to e'/e. In MSSM, the Wilson coef- 
ficients for the AS = 1 effective Lagrangian are modified by integrating out supersymmetric 
particles. (Formulae for the SUSY contribution to the Wilson coefficients are given in Ap- 
pendix 0.) With the SUSY contribution to the Wilson Coefficients given at the electroweak 
scale, we can calculate e'/e using the same prescription as the SM case. 

One important observation is that, in MSSM, AI = | contribution can be largely en- 



hanced relative to the standard model case 0|. This is because the AI = | contribution 



is, as mentioned before, of 0{a^) in the standard model while it can be of 0(ag) in MSSM 
if the mass matrix of the down- type squarks has non- vanishing 1-2 element. Combining the 
fact that AI = | contribution is enhanced by the factor u^'^ ~ 22, the SUSY contribution to 
the e'/e parameter can be order of magnitude larger than the experimental value. Thus, in 
some case, constraint on the CP and fiavor violating parameters in MSSM from e'/e becomes 
severer than that from e. 

As we will see below, the SUSY contribution to e'/e is relatively small in our frame- 
work. Before showing the numerical results, let us explain how the the smallness of e'/e is 
understood in our framework compared to the result given in Ref. [0]. 

As a first step, let us briefiy review how the large SUSY contribution is realized in Ref. 



3^ . The most important enhancement is for the the Li?-type A J = | amplitude by the 



factor of ~ Og/o;^ relative to the SM contributions. The same enhancement also exists for 
the LL- and RR-tjpe operators. However, in taking the relevant matrix elements of the 
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LR-tjpe operators, there is a chirality-enhancement factor proportional to m\/m?g ~ 10, 
and hence the Li?-type operators are more important than the LL- and RR ones. 

Since the A/ = | amphtude is an isospin-breaking effect, a mass difference between 
the up- and down-squark masses is necessary. For the right-handed up- and down-squarks, 
significant mass sphtting is reahzed if the SUSY breaking mass parameters for them take 
different values. On the contrary, the left-handed ones are both from the SU(2)i-doublet. 
As a result, their mass splitting is only from the electroweak symmetry breaking effect 
(i.e., vacuum expectation values of the Higgs bosons), and the mass splitting between left- 
handed squarks is too small to generate significant contribution to the isospin-breaking 
amplitude. Thus, in order to generate Li?-type operators, sizable 1-2 element of the left- 
handed down-type squark mass matrix is necessary. In summary, the condition for the large 
SUSY contribution to e'/ e is (i) mass splitting of between u^i and rf/j, and (ii) non- vanishing 
imaginary part of [m^ ]i2. Combining all of these effects, the SUSY contribution to e'/e can 
be as large as 10~^ when Im([mj ]i2/''^|) ~ 10~^. 

Now, we turn to our case. In our framework, large off-diagonal elements are generated 
for the right-handed down-type squark mass matrix since they are related to the neutrino 
Yukawa interactions. Since sizable mass splitting between up- and down-type squarks is 
possible only for the right-handed sector, only the i?i?-type contribution is modified for the 
A/ = I process, and hence the chirality enhancement is not significant for AJ = | process. 
In addition, if we take the universal boundary condition for the squark masses, mass splitting 
between ur and du is fairly small. As a result, correction to the A/ = | amplitude becomes 
small. 

To make more quantitative discussion, we calculate the SUSY contributions to the AJ = | 



and A/ = | part of e'/e. In our analysis, we followed the prescription given in Ref. I^TI to 



calculate the SUSY contribution to e'/e. We first calculate the SUSY contribution to the 
Wilson coefficients at the electroweak scale. Then, we run them down to the charm quark 
mass scale using the relevant RGEs and took the matrix elements. The resultant SUSY 
contribution is linear in the Wilson coefficients, and numerically we found 

(^y^)SUSY ^ ^ ^^i^ [Cif = mw) + C^Rt'\^^ = mw)] , (3.31) 

i 

where Ci\^~^^ and C]^f~^^ are SUSY contribution to the corresponding Wilson coefficients, 
and the numerical values for Ni are given in Table |I[ 

First, we present the SUSY contribution to the AJ = | and AJ = | part of e'/e normal- 
ized by the standard model ones: 

^,//,nSUSY (j/ASVSY 

^1/2 = 7T7:tsm — ' ^3/2 = TjT-TsM — • y'^-'^'^) 



(^V<=)a/=1/2 (^V^)a/=3/2 



In Figs. § and H, we plot contours of constant R1/2 and R3/2 on rriq vs. (f i — ^2 plane for 
the case of the large angle MSW. Here, we take M,y^ = 3 x 10^^ GeV, and tan/3 = 3. As one 
can see, the SUSY contribution is at most a few % of the standard model contribution for 
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Figure 4: Contours of constant R1/2 in units of 10 ^, which is (eV^)A/l'i/2 normahzed by 
the standard model contribution. Here, we consider the large angle MSW case and take 

, = 0.13 GeV. 



M^^ = 3 X 10^* GeV, tan /3 = 3, and m 



t:^ 0.5 



CM 




600 800 1000 1200 
mq[GeV] 



1400 



Figure 5: Contours of constant in units of 10 ^, which is (eV^)A/l=3/2 normalized by 
the standard model contribution. Here, we consider the large angle MSW case and take 
M^^ = 3 X 10^^ GeV, tan /3 = 3, and = 0.13 GeV. 
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0.11 GeV 


0.13 GeV 


0.15 GeV 


N3 


6.64 X 10^ 


8.92 X 10^ 


1.42 X 10^ 




-9.98 X 10^ 


-9.01 X 10^ 


-8.39 X 10^ 


N, 


1.44 X 10^ 


1.07 X 10^ 


8.29 X 10^ 


N, 


4.30 X 10^ 


3.22 X 10^ 


2.52 X 10^ 


N7 


8.28 X 10*^ 


5.97 X 10^ 


4.49 X 10^ 


Ns 


2.60 X 10^ 


1.89 X 10^ 


1.43 X 10^ 


Ng 


9.08 X 10^ 


1.10 X 10^ 


1.23 X 10^ 


iVio 


2.32 X 10^ 


3.46 X 10*^ 


4.19 X 10^ 



Table 1: Coefficients Ni in units of GeV^ for the ffiting formula of the e'/e parameter. We 
use the bag parameters given in Ref. |pl[, and = 0.11, 0.13, and 0.15 GeV. 



the A/ = i part, and ~ 10 % for the AI = I one. It is notable that, in our simple analysis, 
universal scalar mass is assumed as a boundary condition. Consequently, mass splitting 
between the right-handed up- and down-type squarks becomes small since the masses of 
these squarks are mostly determined by the boundary condition and the RG effect due to 
the gluino loop below the GUT scale which are universal in this case. As a result, masses of 
the right-handed up- and down-type squarks are quite degenerate, and this fact gives another 
suppression factor for the AI = | contribution. Since the contributions to the AI = | and 
AI = I amplitudes are both small, the SUSY contribution to the e'/e parameter is also at 
most a few % level in our framework. 



4 Conclusions and Discussion 

In this paper, we discussed CP violations in the supersymmetric SU(5) model with right- 
handed neutrinos. In this class of model, off-diagonal elements of the right-handed down-type 
squarks are generated via the RG effect even if they vanish at the cut-off scale (i.e., in our 
case, the reduced Planck scale). In general, such off-diagonal elements contain CP violating 
phases and they can be extra sources of the CP violations in the low-energy processes. 
In particular, it was emphasized that such phases can be related to phases in the unified 
theories, which are not related to the parameters in the standard model. 

We paid particular attentions to the CP violations in the kaon system. Most importantly, 
we have seen that the e parameter can be severely affected by the SUSY contribution; even 
with a relatively heavy squark mass of ~ 1 TeV, the SUSY contribution to e can be as large 
as the experimentally measured value if the Majorana mass for the right-handed neutrinos 
is as large as ~ lO^'' GeV. Of course, the SUSY contribution to e strongly depends on the 
phases in the off-diagonal elements of the squark mass matrix. As we have seen, the phases 
in these parameters can be naturally large due to the phases in the GUT model. 
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We have also calculated the SUSY contribution to Br^Ki — > tx^vv) and e'/e. Unfor- 
tunately, however, the SUSY contributions to these quantities are relatively small. In our 
framework, SUSY contribution to Br{KL —>■ n^uu) is a few %, and the SUSY contribution 
to e'/e is also a few % level of the standard model one. 

Thus, in our framework, the SUSY contribution is the most important for the e parameter. 
This fact has an important implication to the future test of the unitarity of the KM matrix, 
since e provides one of the important information of the magnitude of the CP violation in the 
KM matrix in the standard model. In the standard model, the so-called p vs. rj plane was 
constrained so that the standard-model prediction of the e parameter agrees with observed 
one. If the SUSY contribution to e is sizable, however, e cannot be used to constrain the p vs. 
rj plane, which may change the bounds on p and rj. We have seen that the upper and lower 
bounds on the rj parameter changes by about 10 % and 35 %, respectively, if we discard the 
constraint from t]. The deviation from the standard-model prediction on rj may be tested by 
future experiments. In particular, the measurements of 0i and Br{KL — >■ 'k^uv) will provide 
an interesting test of the value of rj. Notice that there are still sizable uncertainties in the 
theoretical calculation of the e parameter, in particular from the bag parameter in taking the 
hadronic matrix elements and from the Wolfenstein's A-parameter in the KM matrix. Thus, 
reduction of the uncertainties in these quantities will be very important to find a signal of 
the SUSY loop using the e parameter. 

Note Added: In finalizing this paper, we found a paper by S. Back et al. |^ which 



has some overlap with our analysis. In particular, in |3^, authors paid special attention to 
non-minimal contribution to the Yukawa matrices at the GUT scale in order to realize a 
realistic unification of the down-type and charged-lepton Yukawa matrices. 

Acknowledgment: One of the authors (TM) would like to thank H. Murayama for useful 
conversations. The work of Y.K. is supported by the Japan Society for the Promotion of 
Science. The work of T.M. is supported by the Grant-in-aid from the Ministry of Education, 
Culture, Sports, Science and Technology, Japan, No. 12047201. 



A Interaction Lagrangian 

In this section, we show the interaction Lagrangian used in this paper. Relevant terms in our 
calculation are the vertices for charginos, neutralinos and gluinos with squarks and quarks. 
At first, we begin relation between gauge and mass eigenstates and mixing matrix, which 
make mass matrix diagonalized. 

The mass and gauge eigenstates for squarks are denoted with (Ja {A = 1, ■ ■ ■ , 6), qiL/pi 
{i = 1,2,3, label for generation), respectively. Then the relations between gauge and mass 
eigenstates are given by: 

QiL = [U{q)]iAqA, qiR = [Uiq)](i+3)AqA, (A.l) 

where U{q) is a 6 x 6 unitary matrix which diagonalize squark mass matrix Ai^, 

[UHq)Mp{q)]AB = mljAB. (A.2) 
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Relation between gauge and mass eigenstate for charginos and neutralinos are given by 



H 



IL 



w 

H. 



R 



2R 

Bl 

wl 

Hi, 



[U{XL)]iXXxL, 
[UiXR)]iXXxR, 

mx%xx'xL, 



(A.3) 
(A.4) 

(A.5) 



where U {Xl/r) ^ (x°) 2x2 and 4x4 unitary matrices which diagonahze mass matrix 
for charginos and neutrahnos, respectively Diagonalization of mass matrices gives masses 



of charginos m-- and neutralinos m^o , 



[U\xr)M^-U{x-l)]xy = m--5xY, [U\x')M^oU{x')]xy = m^ojxv, (A.6) 



For quarks, we chose the basis in which the mass eigenstates {u, d) and gauge eigenstates 
{u', d!) are related as d'j^ — d,, u'^^ — Vkm ul, u'j^ — ur and d'j^ — dR. 

The interaction Lagrangian for chargino-quark-squark couplings are given by 



L 



X-du - Xx [^dAX^L + C^^xPr) 

W{-C^) (C^axPl + C^axPr) Xxd\ + h.c, 



(A.7) 



where C is the charge conjugation matrix, Pl/r = | (1 T 75), and u and d denote up-type 
{u, c, t) and down-type (d, s, h) quarks, respectively. are given by 



^dAX 



r<R 

^dAX 



^uAX 



r<R 

^uAX 



-g2[U*{xR)\ix[U*{u)\aA + [U*{xr)\2x {u)\^^+UyuV^MU (A.8) 

-[U\x-L)]2x[U\u)U[YD]d, (A.9) 
-92[U{xL)]ixY.\U*{d)],A[V^y,U 

-[U{xl)]2xY^!^*^d)]{z+j)A [YdVI^],^, (A. 10) 

[U{XR)]2xjZ[U\d)],A[VlMY^jU (A.ll) 



where §2 is the gauge coupling constant of SU(2) and u,d = (1, 2, 3) for {u, c, t) and {d, s, b), 
respectively. 
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Neutralino-quark-squark coupling couplings are given by 



where 



dAX 



+Xx {KaxPl + N^axPr) uu\ + h.c, 



"^^^'[U{x')Vx + ^[Um2x] [U*{d)U 



[YdU, 



3 



[U*ix')]ix[U*ma-,3)A + [U*ix')hx[U*{d)UA[Y^U, 



uAX 



[U{x%x 



-[U{x%x][U*{u)]uA 



(A.12) 



(A.13) 
(A.14) 



R 

uAX 



-[U{x')UxT.iU*mij+3)A[YuVKMhu, 

^r(x°)]ix[t/*(H)](.+3M 



(A.15) 



-[u*{x')UxY.[u*iu)]jA[viuyM^u, 

with gi being the gauge couphng constant of U(l)y. 

The interaction Lagrangian for gluino-quark-squark coupling is given by 



where 



+u\6-T^ {G^^Pl + G^^Pn) u + h.c, 



^dA 



-^/2g,[U*{d)]iA, 

dA - ^9s[U*{d)\ii+^)A, 



G 



uA 



uA 



-V2g,Y.[U*{u)U[Viu]ju, 

V2gs[U*{u)](^u+3)A: 



(A.16) 



(A.17) 

(A.18) 
(A.19) 

(A.20) 
(A.21) 



with gs being the gauge coupling constant of SU(3)c- Here we take soft-breaking mass 
parameter of gauginos as real positive. 
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B SUSY Contribution to Kl n^iyfy 



In this appendix, we show the SUSY contribution to the Wilson coefficient of effective 
operators for Kl n^vi'. 

There are sets of SUSY loop diagrams in which charged- Higgs, charginos and neutralinos 
are exchanged, which induce effective operators 



Q 



LL 



RL 



(B.l) 



The Wilson coefficients of corresponding operators are referred to C^l and C^^, respectively. 
The effective Hamiltonian is given by Tieff = (Cff'" + Cfr'')QLL + ' + CIT'IQrl- 

In the following, each SUSY contribution of charged- Higgs, charginos and neutralinos are 
separately shown for Z-penguin diagrams. 



LL I 



c 



LL I 



C 



LL 



2Ml 



(B.2) 



dz /^L 



xhxY 



3 



Y.[U*{u)U[U{u)U - Sab (l + 2T|J 



rl/ 2 2 2 \ 



+Sab 



2m,-m^-r(xZ)]ix[t/(xZ)]iy/3K-'<-'^L) 



Xx ~X.Y 



8M| '^^^^^^ 



(B.3) 



Tlf 2 2 2 \ 

/3(m-o^,m^-^,m^-J 



+8ab 



Xx Xy L 6\ ) dA' ^ Xx ' Xy ' 



(B.4) 



with gz = Jgi+ 92 and 



ud 



ttR 
^ud 

or 
or 



g2 cot /? 
g2 tan 



[VKM]„d : 



r(x°)]3x[t/(x°)]3y 



r(x°)]4x[^(x°)] 



4y, 



(B.5) 

(B.6) 

(B.7) 
(B.8) 
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Here T^^ = Tj^ = 0, T^^ = — Tj^ = 1/2. C^^'" is obtained by interchanging R with L, 
The contributions from box-diagrams are given by 



C 
C 



LL I 

Box,!' I 
RL I 



X 



C 



LL 



'^RL I 



-l^x-.'^x-.C^lxC^AYC!;*BxCi;BYl'Arn'^^^^ (B.9) 



2 Xx Xy-a^^ -s^-i- -fs^ -i^-D-i- -'i\- "Xx' ""Xy' 



Xx Xy 



1 ""UA^ Ib'^ 



^^^uBX^^uBY 



Xx Xy 

+NdlYNsAxll i^ll , ' ' <B ) 
1 

4 



N^bxN!:by [2A^f;yiV.V^x5,^x?.^4K^, ^|o., '^L, ^L) 



Xx ' Xy ' " "MA ' > 



(B.IO) 



(B.ii; 



(B.12) 



where C^^x^ ^uAx vertices for x — u — l,x^~^~^^ respectively, that are anafogue to 
those in the quark sector. 



C SUSY Contribution to the AS* = 1 Wilson Coeffi- 
cients 

In this Appendix, we present the formulae for the SUSY contribution to the Wilson coeffi- 
cients of AS" = 1 operators which are used for the calculation of e'/e. To make the notation 
simpler, we first give the formulae for the Wilson coefficients for the following operator basis: 

q 

Hsal.PLdp) E [ctI\qpl,PLqo) + Ctl\qp^,PRq^) 

+ {L^R). (C.l) 

The Wilson coefficients in the basis given in Eq. ( |C.1| ) is converted to the Wilson coefficients 
for the operators Qil^^~^^ — Q[q^~^^ as 

c[T" = ^ti> + lct!>. (C.2) 

CiT" = l-^ClP + lctr. (C.3) 

^(A5=l) _ 1 ^«(S) l^d(S) .p.N 
^L,5 — 12 LR + Q^-R ' V^-V 

^(A5=l) _ 1 ^«(T) l^d(T) P-x 
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MAS= 


=1) = 


1 

J-^«{S) 


1 


r<(As= 

'-'1,8 


=1) ^ 


1^«{T) 




MAS= 
'-'1,9 


=1) = 


1^«{S) 
g^LL - 


gC^LL , 


r<(As= 

'-'1,10 


=1) = 







(C.6) 
(C.7) 
(C.8) 
(C.9) 



and 4r = 



(AS=1)| 



Let us first consider the effect of gluon-penguin operator wliicli contributes only to the 
A/ = i amplitude. Denoting 



12N. 



j2( 2 2 2 2 2\ 



c t2/ 2 2 2 2 2 \ I -"c 2 7-1/ 2 2 2 2 2 \ 



6 

- ^dsC^lxCdAX^^ , ^Ia ' ™Ia ' "^L ' ™Ia ) 

1 TirL* TirL T2f 2 2 2 2 2 \ 

-Q93N^AxNdAxh (^xl^^dA^^dA^^dA^^dA^^ 



we obtain 



C 



QiS) 
LL 



GP 



"?(S) 



1 ^GP ^q(T) 



'LL 



GP 



= 93CY 

GP 2^ ^ 



GP 



(C.IO) 



(C.ll) 



Notice that, in Eq. (|C.10|) and hereafter, summation over the dummy indices [A, X, and so 
on) is implied. 

Contribution of the photon-penguin diagrams is parameterized by 



-1 

3gj\r dA 5\ G' dA' 0!^' ^a' ^A' 



-'-r2/2 2 2 2 2\ 

-L[m-±,m- ,m- ,m-- ,m-- 

n 5 V ' MA ' MA ' "A ' "A 



I 72/2 2 2 2 2 



Xx' x^' xj' x^' 
"2rl/2 2 2 2 2\ 

--m-± it- (m-± , m-± , m-± , m-± , m,- ^ ) 



3 
1 
2 
1 



+ Y^eiVi;,*^iV,V/5^(m|o^, m|^, m|^, m|^, m|J. 
With the photon-penguin contributions to the Wilson coefficients are given by 



C 



■9(S) 
LL 



PP 



c 



g(s) 



PP 



—ee„C 



PP 
q^L ' 



■<?(T) 
LL 



PP 



LR 



PP 



0. 



(C.12) 



(C.13) 
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Finally, we present the box contributions. First, the box contributions with internal 
gluino and/or neutralino lines are given by 



C 



5(S) 
LL 



(G,X°) 
Box 



(G,X°) 
Box 



r<9(s) 



(G,X°) 
Box 



+ 



2 ™2 



71/ 2 2 2 2 \ 

/4(m^,m^,mj^,m,~J 



9s > 



1 



c' 2 7-0/2 2 2 2 

(m^, m^, m^-^, m^-^ 



L* 



1 

4 
1 

2 

L*/^L /^L* /^L 



tL* 



7-1/ 2 2 2 2 \ 

74 (m^,m^,mj^,m^-^) 



(C.14) 



sA^dA^qB^qB 
1 



4iV. 



2 7-0/2 2 2 2\ 



( r^L-^ j^jij r^L^ AT^ _L AT^"^ A/"^* /^-'^ 



~*^sA*^dA*^qB*^qB 



Xx ' <^A ' 



(C.15) 



- 7-1/ 222 2 \ 



2 7-0/2 2 2 2\ 
■g^^G^l^G' "^G' ^d^' 



1 



( /~tL* ]\tL AT^* /^R I /vr^* riH* AT^ 

sA^^ dAX^^ qBX^ qB + ^^sAX'-^dA'-^gS^^gSX 



ri?* 



atR 



I AT-i* AT-i AT-i?* A7-R 7-O / 2 2 2 2 \ 



(C.16) 



25 



c 



</(T) 
LR 



Box 



^ sA^ dA^ qB^ qB 



8N, 



j\i 2 2 2 2 \ 



7V^ — 2 

■"'c ^ 2 7-0/2 2 2 2n 



^-^^ sAX^ dA^ qB^ qBX + ^ sA^ dAX^ qBX^ qs) 

xmam^oJ^{ml, m|o^, m|^, m?^). 
In addition, contributions with internal chargino lines exist, which are given by 



(C.17) 



C 



u{S) 
LL 



^«(S) 



LL 



c 



d(S) 
LR 



(X±) 

Box 

(X±) 

Box 

(X±) 

Box 

(X±) 

Box 



Xx' Xy' dA' 9s > 



-^sBX^dBYGuBX^uBY^l ("^Xx ' "^Xy ' ''"L ' ''"L ) ' 



~:^GsBxGdBYGdBYGdBX^li''^i(t^^^^t'''^lA^''^^ 



CsBxCdBYCdBYCdBxm^^m-±I^{m-±,m-±,ml^,m^ 



Xx ' Xy ' "A ' ■ '"Us ^ 



(C.19) 
(C.20) 
(C.21) 



while Cf7^ and cfP 



Box 



LR 



(X±) 

Box 



vanish. 



The Wilson coefficients and are obtained from and Cl,^ by interchanging 
the indices L and R. 



D Master Integrals 

In this Appendix, we present the master integrals used in the calculations of the loop dia- 
grams. 

The function is defined as 

/22 "2 \ f ^ ^ / 1 \ 

I J, {m„m„ • • • , m^) = y _ ^2)(^2 - rnl) ■ ■ ■ {k^ - ml) ' ^ ' 
Explicit formulae of used in our calculations are as follows: 

4(Mi, Mi, Ml, Mi, mi, mi) = ^ + ^ ^ ~ + ^ " + (D.2) 

48 TT^ Mx (a; — 1) 

Tn/,^2 ,^2 ,^2 ,^2 2 2 n 17 - 9a; - 9a;^ + + 6 (1 + 3a;) logx _ ^, 

7°(Mi, Mi, Mi, Mi, mi, mi) = „^ , , . sr (D-3) 

/|(Mi, mi, mi, mi, mi) = ~ ^'^J , (D-4) 

96 TT^ Mx (x — 1) 
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X 



ll{Ml 
I'M 



2 2 2 2 \ 



2 + 3x-6x^ + + 6x logx 



1 — + 2x log X 
16 7r2M|-(x- 1) 



3 ' 



2 2 2 N 2 - 2x+ (1 + a;) logx 



16 7r2 Mx^ (a: - 1) 



3 ' 



,^2 2\ + x - x logx) 



16 7r2(a;- 1) 
1 — X + log X 
167r2Mx^(;7;- 1)' 

2 , — 1 + X — X logx 
"^a) = ,^ _9 /„ ^ -'Div, 



167r2 (x- 1) 



2 2 \ 



1 



1 

167r2Mi(x-|/) 
1 



-1 + X — x^ logx 



(x - 1)^ 
— 1 + X — X logx 



— X 



y) 



167r2(x — y) 
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(x - 1)^ 
-X + x^ — x^ log X 



— X 



167r2Mx^ {x-y) 



2 2 2 \ 

mj^,mB,mc) = 
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X — 1 
— X logx 
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(x y) 



y) 

+ -^Div, 



X 



16 tt"^ Mx {x — y) [x — ^ 



167r2M|-(x-y) \x 



-(x^y) 
— x^ logx 

—X logx 



(x ^ 



X z] 



(x 



X 



(D.5 
(D.6 
(D.7 
(D.8 
(D.9 
(D.IO 
(D.ll 
(D.12 
(D.13 
(D.14 

(D.15 
(D.16 



where x = m^/M^, |/ = m%/Mx, z = rriQ/Mx, and Jdiv contains divergence of 1/e in 
dimensional regularization: 



^niv — 



167r2 



1 n Mtt 

+ Ml 



(D.17 
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